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Abstract 

Using modified gravity with non-linear terms of curvature , and 
_R(2+r) ^ being a positive real number and R being the scalar 

curvature), cosmological scenario, beginning at the Planck scale, is ob- 
tained. Here a unified picture of cosmology is obtained from f{R) — 
gravity. In this scenario, universe begins with power-law infiation fol- 
lowed by deceleration and acceleration in the late universe as well as 
possible collapse of the unverse in future. It is different from f{R) — dark 
energy models with non-linear curvature terms assumed as dark energy. 
Here, dark energy terms are induced by linear as well as non-linear 
terms of curvature in Friedmann equation being derived from modified 
gravity . It is also interesting to see that, in this model, dark radiation 
and dark matter terms emerge spotaneously from the gravitational sec- 
tor. It is found that dark energy, obtained here, behaves as quintessence 
in the early universe and phantom in the late universe. Moreover, anal- 
ogous to brane-tension in brane-gravity inspired Friedmann equation, 
a tension term A arises here being called as cosmic tension. It is found 
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that, in the late universe, Friedmann equation (obtained here) contains 
a term — p^/2A (p being the phantom energy density) analogous to a 
similar term in Friedmann equation with loop quantum effects, if A > 
and brane- gravity correction when A < 0. 



1. Introduction 



Cosmology was revolutionized by observations made during last few 
years [1, 2]. These observations show conclusive evidence for acceler- 
ation in the late universe, which is still a challenge for cosmologists. 
Theoretically, it is found that dark energy (DE) violating strong en- 
ergy condition{SEG) or weak energy conditionCWEC) is responsible for 
it. So, in the recent past, many DE models were proposed to explain 
the late cosmic acceleration. A comprehensive review of these mod- 
els is available in [3]. Later on, it was realized that even non-linear 
terms of curvature with n > also could be used as DE [4]. 
Although this model explained late cosmic acceleration, it exhibited 
instability and failed to satisfy solar system constraints. It was im- 
proved further by Nojiri and Odintsov taking different forms of f{R) 
for DE. These improved models satisfied solar system constraints ex- 
hibiting late cosmic acceleration for small curvature and early inflation 
for large curvature. Thus, in f{R)— dark energy models, non- linear 
curvature terms are considered as an alternative for DE [5, for detailed 
review]. Recently, in [6], it is shown that f{R)— dark energy models 
with dominating powers of R for large or small R can not yield viable 
cosmolgy as results contradict the standard model and do not sat- 
isfy Wilkinson Microwave Anisotropy Probe (WMAP) results, though 
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these models pass solar system constraints and explain late accelera- 
tion. Moreover, it is also shown that the most popular model with 
f{R) ^ R + aK^ + PR~'^{m > 0, n > 0) considered in [5] is unable 
to produce matter in the late universe prior to the beginning of late 
acceleration [6]. 

In what follows, it is aimed to get a viable cosmolgy consistent with 
WMAP from f{R)— gravity, where f{R) = R/I671G+ powers of R, not 
from f{R)— dark energy models [7, 8, 9] discussed above . There is a 
crucial difference between the two. In the latter case, which is criticised 
in [6], non- linear terms of curvature are treated as dark energy terms. 
On the contrary, in the former case, neither linear nor non-linear term 
is considered as dark energy. In the present model, it is important to 
see that DE terms are induced by linear (Einstein-Hilbert term) as well 
as non-linear terms R^ and R^^'^'^^ in the action. In f{R)— dark energy 
models, dark energy terms depend on f{R) terms and its derivative 
F = df /dR. In the former case, induced DE terms depend on the scale 
factor a{t) of the homogeneous and flat of Friedmann - Robertson - 
Walker (FRW) model of the universe. 

In this paper, the f{R) — gravity based modified Friedmann equa- 
tions are derived in the early and late universe taking small and large 
a{t) respectively. Contrary to f{R)— dark energy moofe/s dominance of 
non-linear terms of curvature is not taken here for small and large R as, 
in the present model, DE terms emerge as imprints of both linear term 
as well as powers of R. In f{R)— dark energy models, action contains 
lagrangians for matter and radiation [5, 6]. It is interesting to see that, 
in the present f{R) — gravity model of cosmology, dark radiation and 
dark matter terms emerge spontaneously [8, 9]. 
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Friedmann equation (FE), obtained here, gives cosmic dynamics, 
where some terms emerge having forms different from known forms 
of energy density (radiation or matter) and violate SEC in the case 
of early universe and WEC in the case of late universe. So, these 
terms are recognized as curvature induced DE. Thus, this approach of 
getting DE from curvature is different from the approach of [5, 6]. It 
is also interesting to see that, in the case of late universe, FE contains 
square of DE density with negative sign analogous to a similar term in 
Friedmann equation with loop quantum effects [10]. 

In what follows, it is interesting to see that radiation and matter 
terms in FE (obtained here) emerge spontancosly if r = 3 and n = 
1/4. In [5], radiation and and matter terms do not emerge from the 
gravitational sector. It is important to mention that, here, theory 
itself ignores the cases of r = 1, 2 i.e. theory suggests that non-hnear 
of R should be i?^ and R^. In [5], a certain form of scale factor a{t) 
is assumed and, later on, conditions are obtained for assumed a(t) 
consistent with experiments and satisfying stability criterion. Here, 
a{t) is not assumed, but it is derived solving Friedmann equations in 
the different stages of the universe. This is another important difference 
in approach of this paper compared to [5] . All these results are obtained 
from modified gravity without using any exotic matter or field. This 
approach is adapted in [7, 8, 9] also. 

In [9] , a unified picture of the universe, from early inflationary stage 
to late acceleration and deceleration driven by dark radiation and dark 
matter between these two stages, is obtained taking linear and non- 
linear terms of curvature as well as a different scalar. But, in this 
paper, the same result is obtained from curvature terms only . 
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Investigations, given below, show that, originating at the Planck 
scale, the universe inflates for a very short period followed by decel- 
eration driven by curvature-induced dark radiation and subsequently 
by dark matter. In the very late universe, around 12.86Gyr, a transi- 
tion from deceleration to acceleration takes place. Further, it is found 
that late acceleration will continue upto 58.48Gyr. This is an epoch for 
transition from acceleration to deceleration. At this epoch, acceleration 
will stop and deceleration, driven by matter, will resume [9, 12]. The 
decelerated expansion will continue upto the time ~ 1.05 x lO^^^Gyr. 
By this time, the universe will have maximum expansion. So, it is 
natural to think that the universe will retrace back and contract. Re- 
sults show that, due to contraction, universe will collapse by the time 
~ 1.3 X lO^^^Gyr. 

In [13], quintessence DE in the early universe and phantom DE in 
the late universe have been considered taking non-gravitational scalar 
field (curvature independent scalar field) as DE source. The present 
paper is different from [13] in the sense that here we have gravitational 
origin of quintessence and phantom DE as it is obtained in references 



Natural units {ks — h — c — 1) (where kB,h,c have their usual 



6.58 X W-^^sec. and IGeV = 1.16 x 10^^^. 



2. Action for f{R) — gravity and Friedmann equations 



Here action is taken as 



[7, 8, 9]. 



meaning. GeV is used as a fundamental unit and we have IGeV ^ = 




(2.1) 
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where G = Mp'^{Mp = lO^^GeV is the Planck mass), a is a dimension- 
less coupling constant, P is a, constant having dimension (mass) (as 
R has mass dimension 2) with r being a positive real number. 

Using the condition SS/Sg^'' — 0, action (2.1) yields field equations 

1 1 „ . _„ 1 ^2 



+/3(2 + r)(V^V.i?(^+'-) - g,.nR^'^^)) + ]^Pg,,R^^+^^ 

-P{2 + r)R^'+'^R^, = 0, (2.2) 

where V^^ stands for the covariant derivative. 
Taking trace of (2.2), it is obtained that 

^ - GaOR -3p{2 + r) + prR^^^''^ = (2.3) 



IGttG 

with 
In (2.3) 



□^(i+r) _ (i + r)[i?"ni? + ri?('^-^)v^i?V^i?]. (2.5) 



Prom (2.3) and (2.5) 
R 



^ ■[6a + 3/3(l + r)(2 + r)i?']n7?-3/3r(l+r)(2+r)i?("-^)v'^i?V^i? 

+/3ri?('+'') = (2.6) 

In (2.6), [6a + 3/3(1 + r) (2 + r)i?''] emerges as a coefficient of OR due 
to presence of terms aR^ and jSR^"^^'^'^ in the action (2.1). If a = 0, 
effect of R'^ vanishes and effect of R^'^^^^ is switched off for (3 — So, 
like [9] an effective scalar curvature R is defined as 

7i?^^- = [6a + 3/5(1 + r) (2 + r)i?n, (2.7) 

where 7 is a constant having dimension (mass)"^'' being used for di- 
mensional correction. 
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Using (2.7) in (2.6), we have 

_ (7/r)^^y(V--2)[yny + (l/r - i)v'^y v„y] 

IdttG 

-3(/3/r)(l + r)(2 + r)y(^/'^-^)v'^yv^y + /3ry(2+'')/'' = 0, (2.8) 



where 

(2.8) is simphfied as 
1 



7i?^ - 6a 



3/3(1 + r) (2 + r)' 



(2.9a) 



Y - {-f/r)ErpY + (l/r - i)y-^v'^yv^y] 



-3(/3/r)(l + r)(2 + r)V''yv^y + /3ry(^/^+^^ = 0. 



(2.96) 



Using (2.9a)in (2.9b), it is obtained that 



167rG'7 



6a 



-M-(l/r-l) 



7 



+[3/3Ml + 0(2 + 0/7']^'' 
(2.10) is re-writtcn as 

6a 



[6a - -fR"' 

-fR' - 6a 



■y^'Rv aR+R'''y^Ry uR" 



3/3(1 + r) (2 + r) 



(l/r+2) 



= 0. (2.10) 



1671^ ^R-^ 



-il-r) 



+ nR + {r- 1)^-^V^^V^^ 



^R 



r-l 



6a — jR 



-v^Rv^R + rR-^y^^Ry^R 



+r)(2 + r)/Y]R 



,21 52r--l 



7i?'' - 6a 



(l/r+2) 



= 0. 



(2.11) 



.3/3(1 + r)(2 + r) 

Experimental evidences [14] support spatially homogeneous flat model 
of the universe 



dS^ = dr - a\t)[dx^ + dy^ + dz^ 



(2.12) 



with a{t) being the scale factor. 
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For a{t), being the power-law function of cosmic time, R ~ a"". For 
example, R ~ for matter-dominated model. So, there is no harm 
in taking 

R^—, (2.13) 

where n > is a real number and A is a constant with mass dimension 
2. 

Connecting (2.11) and (2.13), it is obtained that 

2-n-n(r-l)+^ nr - = — 1 

^-1/3 

"n(7>l^a— )'[3r(l + r)(2 + r)] " 7^'■a-'^1'+'/^ (2-14) 
taking — —(3'^/^ and ignoring complex roots as these roots 

lead to unphysical situations. Now, we have foUwing two cases. 

Case 1 : The Early Universe 

In this case, a{t) is very small, so (2.14) is approximated as 



a 

-+ 
a 



- + 
a 



2 — n — nr 



a) 



n(7A'-a-"0^[3r(l + r)(2 + r)]i+V'- 
x[6Q;-7A''a-"n'^'^' (2-15) 



as 

7A'a" 



-1. (2.16) 



6q; - 7A'-a-'*^ 
Integration of (2.16) leads to 

/d\2_ B 2(3-^^1'- 

\a) ~ a(2+2A^) ~ n{-iA'-f[?>r{l + r)(2 + r)]i+i/r-a(2+2M) 

X y"a(™+2"'-)[6a-7A^a-"'-]2+Vr- (2.17) 

with 

M = 2-n-nr. (2.18) 
Case 2 : The Late Universe 
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In this case, a{t) is large, so (2.14) is approximated as 



+ 



2 — n — n{r — 1) — nr ^-^ ~ 



- 1 



n(7A'a-"'-)2[3r(l + r)(2 + r)]i+iA 



(6Q;)'+^/'-[a2"" - (2 + l/r)7A'^a"''] 

(2.19a) 



as 







6q; — 7 A'' a 

for large scale factor a. So, (2.19a) is re- written as 



a 

- + 
a 



2 - 2nr 



d\ ■ 
a) 



(2.196) 



where 



/ 6a X r 1 _ [3r(l + r)(2 + r)]-i-V- / 6a xi 

^" l7A-ni67rGn + ^ ^ l^^^.j 



7^^ 



(2.20a) 



and 



E 



/ 6a \2 



.7A 



1 [3r(l + r)(2 + r)]-i-V'- / 6 



leTrGn 



(2.19b) is integrated to 



(-) =- 

\a/ a 
with 



B 



2D 



(2+2JV) (2 + 2A^ + nr) L D{2 + 2N + 2nr) 



E{2 + 2N + nr) 



(2.21) 



N = 2- 2nr. (2.22) 

Further, it is found that if M = 1, the first term on r.h.s. (right hand 
side) of (2.17)gives radiation. Moreover, UN — 1/2 the first term of 
r.h.s. of (2.21) gives matter. So, using M = 1 in (2.18) and N = 1/2 
in (2.22), it is obtained that 

3 
4' 
1 



nr — 



n 



(2.23) 
(2.24) 
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and 
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(2.25) 



3. Power-law inflation followed by deceleration 
in the early universe 

The approximated Friedmann equation (2.17), in the case of the 

early universe, looks like 

2 B 8/3-^/3 



/ay _ B 



z9/2[6Q;-7^=^a-=^/T/' (3.1) 



(7^3)2 [l80]4/3a4 

using definitions of M and N as well as (2.24) and (2.25). In (3.1), 



j a9/2[6a-7^3a-3/4]7/3 



11 



3^-3/417/3 



7 

22 



X j a}''l\Qa - 7/l3a-3/4}^/3^a. (3.2a) 



It is noted that for 



\ 4/3 



a < (^-fA^/6aj ' = a^, (3.26) 



terms within bracket and the integral on the right hand side of (3.2a) 
are of the order of a^^/^. 
So, 



/ 



a9/2[6a - 7A3^-3/4]7/3 ^ [AalV2|6a _ 7^3^-3/4|7/3_ 



Thus, using (3.2a,b,c),(3.1) is approximated as 

2 B IQp-^/^ 



3/2 



a-3/4 _ ^-3/4 



7/3 



(3.2c) 



(3.3) 



ll(7A3)-l/3[180]4/3 

It is interesting to sec that a radiation density term B/a'^ emerges 
spontaneously. This type of a term, being called dark radiation, emerges 
in brane-gravity inspired Friedmann equation too. So, analogous to 
brane-gravity, here also B/a'^ is called dark radiation. Other terms on 
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r.h.s. of (3.3) arc caused by linear as well as non-linear terms of cur- 
vature in the action (2.1). These terms also constitute energy density 
term 



qu 



16/5-'/' 



7/3 



87rGLll(7A3)-V3[i80]4/3 JL" ^A^\ ^ ' 

(taking real root of (—1)"^/^ as above) satisfying the conservation equa- 
tion 



Pde + 3-(pde + Pde) = 0. 



(3.5; 



Connecting (3.4) and (3.5), equation of state (EOP is obtained as 



where 



/ 



I6/3- 



-1/3 



(3.6a) 



(3.66) 



87rG'll(7A3)-i/3[i8o]4/3 
(3.6a) is the scale factor-dependent equation of state parameter, valid 
for Op < a{t) < ttc- Such an equation of state parameter is obtained in 
[8] also. It yields 

Pde + Pde > o and pf^ + 3p^" < o, 

for ap < a{t) < a^. It shows that DE, having energy density (3.4) 
mimics quintessence dark energy [7, 8, 9]. 

Here investigations start at the Planck scale, where DE density is 
obtained around lO^^GeV^. So, (3.4) is obtained as 



-3/4] 



7/3 



with 



75^-3/2 



10'°a 



-3/4 _ ^-3/4 



-7/3 



(3.7a) 



(3.76) 



Thus , (3.6b) and (3.7b) imply 



(3.7c) 
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Connecting (3.3) and (3.7a), it is obtained that 



ay 
a) 



B 



Stt X lO^V a \3/2ra-3/4 



ap) 



-3/4 



-3/4 



— a< 



-3/4 



7/3 



(3.8) 



using G = = 10-38GeV-2. 

It means that the universe is driven by radiation for a > Uc- More- 
over, (3.7a) shows that vanishes at a = Oc and for op < a{t) < Oc, 
cosmic dynamics is given by 



Stt X 10 



37 



a \3/2 

ap) 



-3/4 



-3/4 



-3/4 



-3/4 



7/3 



Stt X 10 



37 



a \ -1/4 

Op/ 



(3.9) integrates to 

a{t) — ap 



1 + 10 



18 . 



— it-tp) 
12^ ' 



(3.9) 



(3.10) 



showing acceleration as a > 0. 

If expansion (3.10) yields sufficient inflation in the early universe. 



^ = 10^^ 

ap 



(3.11) 



The universe comes out of the inflationary phase at t = tc when a{t) 
acquires the value Oc. So, from (3.10) and (3.11), it is obtained that 



tc^tp + m 



-18 




ap) 



2.76 X lO^tp 



(3.12) 



using (3.11). 

For a > ac, we have Priedmann equation (3.8) as 

(-)' - 

\a/ a^ 

This equation integrates to 

a{t) = ac[l + VB{t - QY^^. 



B 



(3.13) 



(3.14) 



(3.14) yields a < showing deceleration driven by dark radiation term. 
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4. Deceleration followed by acceleration in the late universe 
as well as future collapse of the universe 

In the late universe, the effective Priedmann equation is given by 
(2.21). Using (2.23)-(2.25) in (2.21), we obtain 



ay 
aJ 



C 8D 



+ 



,3/4 



15 



5E 



where 



and 



-3/4 



E 



AttG 

-3/2 

ac 



287rG / 1 
135 U80/3 



6^" 

1/3 



3/4 



a. 



-1/4 



1 - 



AttG / 1 \V3 



135 V 180/3/ 

being obtained from (2.20a) and (2.20b) using (2.23)-(2.25) and (3.2b). 

The first term, on r.h.s. of (4.1a), emerges spontaneously and has 
the form of matter density, so it is recognized as dark matter density 
like dark radiation. Moreover, the second and third terms on r.h.s. of 
(4.1a) emerges due to linear and non- linear terms of curvature. It is 
interesting to see that if 



-1/4 



(4.1a) 



(4.16) 



(4.1c) 



3/4 



and 



(4.1a) looks like 



a) 



D 

5^" 

3D^ 
2bT:GE' 

2 87rGr 3C 



Pde 



A 



87rGa3+^^c{^ 2 A } 



(4.2) 



(4.3) 



(4.4) 



(4.5) 



Conservation equation (3.5) for yields 

ph 5 

w, = . 

dc ^ 

(4.5) shows that the curvature-induced energy density mimes phan- 
tom dark energy as w^g < — 1. Thus, in the late universe, a phantom 
model is obtained from curvature without using any source of exotic 
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matter. Apart from this, (4.4) contains a term —{p^^Y/2\ analogous 
to brane-gravity correction to the Friedmann equation (FE) for nega- 
tive brane-tension [11] and modifications in FE due to loop-quantum 
effects [10]. Here A is called cosmic tension [7, 8, 9]. 

According to WMAP results [15], present density of pressureless dark 
matter is obtained to be Pq"^^ = 0.23^0^^ and present dark energy density 

p£o = o.rsp^o"" with 

^° " SttG' 

where current Hubble's rate of expansion Hq = 100 hkm/M pes econd ~ 
2.32 X 10-^^hGeV and h = 0.68. Thus, 

= 2.9 X 10-^^GeV^ (4.6) 

Using these values, it is obtained that 

3C _ 6.67xl0-^« 



and 



= ^aV4 ^ 2.117 X 10-^^a=^/^ (4.8) 



from (4.2). 

Connecting (4.4), (4.7) and (4.8), it is obtained that 



(°)' = r6.67 X 10-« ^^^^ ^ io-«„3/4 1 1 _ 2.117 X IQ-^a^'*^ 

(4.9) 

(4.9) shows that 

6.67 X 10-^^ 



> 2.117 X 10"^ 
a'' 

for a < 0.735 and 



6.67 X 10-^^ .7 

- < 2.117 X 10-^V/^ 



a3 



for a > 0.735. 
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It means that a transition from matter-dominance to DE-dominance 
takes place at 

= 0.735 (4.10) 

giving red-shift 

= — - 1 = 0.3607 (4.11) 

which is very closed to lower limit of 2;* given by 16 Type supernova 
observations [2]. Thus, for a < 0.735, (4.9) is approximated as 



<a/ 3 
which integrates to 



d\2 StiG [6.67 X 10-^^] 5.59 x 10"^^ 



(4.12) 



a{t) = ad[l + 7.48 x 10-^3^//'(i - td)f^ (4.13) 

It shows decelerated expansion as a < 0. 
When a > 0.735, (4.9) is approximated as 

(3 = = 177 X lO-a^'^ [1 - . (4.14) 

(4.14) integrates to 

[2.117X10-^^ f / _3/4 2.117 X 10-47 ^l21-4/3 

""^^^ = [ 2A +|V"* 5xl0-^3(^_^^)| J _ 

(4.15) 

This scale factor yields a > showing acceleration in the late uni- 
verse. In (4.15), t^ is the time of transition from deceleration to accel- 
eration in the late universe. 

Using the present age of the universe to — 13.7Gyr = 6.6x lO^^GeV"^ 
and Oo = 1 (as given above) in (4.15), is calculated as 

t^^to- 0.4 X lO^^GeV-^ = 6.2 x lO^^GeV"^ = 12.86Gyr. (4.16) 
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(4.14) shows that accelerated expansion (4.15) stops at a = Oe satis- 
fying the condition 

2.117 X 10-^^af/^ = 2A. (4.17) 
a{t), given by (4.15), acquires the value Og by the time 



_3/4 2.117 X 10-47 / _3/4 2.117 X 10-47- 

Oe ^ A/a* - 



/ _3/4 2.117 X 10-47-, 
(4.18) 

The Priedmann equation (4.9) shows that, a,t t > te, expansion of 
the universe is driven by matter again and it reduces to (4.12) yielding 
the solution 

a{t) = ae[l + 1.21 X lQ-^^af^{t - Qf'^. (4.19) 

which shows decelerated expansion. Thus another transition from ac- 
celeration to deceleration will take place ait — te- 
It is interesting to note that the term 

2.117 X IQ-^V/'i- 



2.117 X 10-^ V/^ 



1 - 



2A 

will be negative as a > Cg. So , gradually universe will reach a state , 
where scale factor a{t) acquires its maximum value a^. Ki a — a^, a — 
in (4.9) and satisfies the condition 

6.67 X 10-« ^ ^^^^ ^ 10-«„yf -"^ - l|. (4.20) 

am gives the maximum expansion, so for t > tm universe will change 

its direction and retrace back leading to contraction. As a consequence, 

for t > tm,(i{t) will decrease and a'^ term ,in (4.9), will dominate 

yielding the effective equation 

d _ 8.099 X 10-43 
a a^/^ 
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Here, ^ 



< due to contraction. This equation is integrated to 



a(t) = a^[l - 1.21 X 10-'\-'/\t - tm)]'/\ 



(4.21) 



(4.21) shows that at time 



t = tm + 8.26 X lO^^a^^ 



(4.22) 



a = 0. It means that universe will collapse at this time. 



5. Summary 



Results, obtained above, arc summarized as follows. Here f{R) — 
gravitational action is obtained by adding higher-order terms and 
7^(2+r) Qf scalar curvature R to the Einstein-Hilbert term. Gravitational 
field equations are derived from this action. Using R ~ a~" in trace 
of f{R)— gravity field equation, Priedmann equation is obtained. It is 
found above that if r = 3 and n = 1/4, FE (obtained here) contains 
quintessence like dark energy term as well as radiation like term in the 
early universe. Here, radiation emerges spontaneously and is termed as 
dark radiation which is analogous to a similar term in brane-gravity- 
based FE. Dark energy term is induced by curvature and it vanishes 
when the scale factor a{t) acquires a finite value Oc given by (3.2b). It is 
interesting to see that, in the late universe, dark matter term emerges 
spontaneously and in the very late universe (when it is 12.86Gyrs old) 
the universe is dominated by curvature-induced phantom dark energy 
with w = —1.25. Contribution of R to DE is a physical concept in 
addition to its usual role as a geometrical field. Thus, dual roles of R 
(as a physical field as well as a geometrical field) [16] are manifested 
here. The cosmological scenario, obtained here, from f{R)— gravity 
with higher-order terms R^ and R^, is given as follows. 
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It is found that the early universe inflated for a very short period with 
power-law speeded-up expansion, driven by curvature-induced quintes- 
sence DE. When the scale factor increased upto 10^^ times the scale 
factor at Planck scale, quintessence DE density vanished. As a con- 
sequence, universe came out of the inflationary era. Later on, early 
universe decelerated as t^'''^ driven by dark radiation. Subsequently, 
when the universe became sufliciently old, it decelerated as t"^^^ driven 
by dark matter. At red-shift 2;** = 0.3607, a transition from decelera- 
tion to acceleration took place and universe began to accelerate driven 
by curvature-induced phantom DE explaining the present acceleration 
of the universe. It is found that the late acceleration is transient and 
it stops when phantom DE grows to a flnite value equal to 2 A. Here 
A is the cosmic tension analogous to brane-tension, which is explained 
above. Interestingly, the phantom model (obtained here) is free from 
the menace of future-singularity. Here, it is shown that universe will 
reach its maximum expansion at a time tm- Later on, it will contract 
and collapse in a flnite future time. Thus, it is found that contrary 
to /(i?)— dark energy models, we get a viable cosmology from f{R) — 
gravity. Also, it is found that curvature induced phantom DE has 
a crucial role in the dynamics of present and future universe giving 
different phases mentioned above. Results, obtained above, support 
observations made so far. 
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